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Abstract 

The phase structure of flavoured M = 2 SYM on a three sphere in an external magnetic field 
is studied. The pairing effect of the magnetic field competes with the dissociating effect of the 
Casimir energy, leading to an interesting phase structure of confined and deconfined phases 
separated by a critical curve of a first order quantum phase transition. At vanishing magnetic 
field the phase transition is of a third order. For sufficiently strong magnetic field, the only 
stable phase is the confined phase and magnetic catalysis of chiral symmetry breaking is 
realized. The meson spectra of the theory exhibit Zeeman splitting and level crossing and 
feature a finite jump at the phase transition between the confined and deconfined phases. At 
strong magnetic field the ground state has a massless mode corresponding to the Goldstone 
boson associated with the spontaneously broken U(l) /^-symmetry analogous to the rf meson 
in QCD. 
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1. Introduction 



The influence of magnetic fields on flavoured gauge theories has been extensively studied in 
the literature. In the pioneering works of refs.[T]-[6], it has been shown that magnetic fields 
act as a strong catalyst of mass generation and chiral symmetry breaking leading to the 
formation of a fermionic condensate, even in the slightest attractive potential. The essence 
of this effect is the dimensional reduction, from D to D-2 dimensions, in the dynamics of 
fermion pairing. The effect has been shown to be model independent and therefore has a 
universal nature. Given such universal nature, it is natural to study this effect in holographic 
gauge theories, where one can rely on the powerful techniques of the AdS/CFT correspon- 
dence at strong coupling. 

In its original formulation [7], the correspondence relates perturbative string theory on an 
AdSs x S 5 background to a four dimensional M = 4 SYM theory. An important extension 
of the correspondence [8] , making it relevant to the description of flavoured Yang-Mills the- 
ories, was the inclusion of fundamental matter via the introduction of flavour branes in the 
so-called probe limit, where the number of different flavours is much less than the number 
of colors, Nf <C N c . This corresponds to the quenched approximation on the gauge theory 
side and the probe approximation on the supergravity side of the correspondence. 
The holographic approach to magnetic catalysis was initiated in ref. [5], where the holo- 
graphic gauge theory dual to the D3/D7 intersection has been analyzed. Further relevant 
studies have been performed in refs. [10j-[53j. Holographic studies of backreacted flavours 
in external magnetic fields have been performed in refs. [10j-[12j. 

In the present work, we are interested in studying the effect of an external magnetic field on 
a flavoured gauge theory defined on a compact space. A canonical example of holographic 
gauge theory on compact spaces is M = 4 SYM theory on a three-sphere which is dual 
to string theory defined on AdSs x S 5 in global coordinates. The addition of fundamental 
flavours is achieved in the same way as in the flat case, namely by introducing probe D- 
branes. The confinement /deconfinement phase transition of flavours on S 4 has been studied 
in refs. [T3"f 14)115] . and the effect of i?-charge and isospin chemical potentials was addressed 
in refs. [THIfTTj . 

The finite volume gives and extra energy scale associated to the Casimir energy of the theory. 
Despite the fact that our study is at zero temperature, the expected effect of the Casimir 
free energy is to favour the dissociation of meson-like bound states, triggering a confine- 
ment /deconfinement phase transition [14J. On the other hand, magnetic fields favour the 
formation of bound states. Therefore, one would expect to find an interesting phase struc- 
ture of confined and deconfined phases separated by a critical curve across which a quantum 
confinement/deconfinement phase transition takes place. Note that this phase transition 
does not have an analogue in the theory on a flat space. However, for sufficiently strong 
magnetic fields the effect of the Casimir energy will be subdominant and the theory should 
be qualitatively similar to the flat case. Thus we expect that for strong magnetic fields the 
phase transition seizes to exist, and the theory is in a confined phase where the vacuum at 
zero bare mass spontaneously breaks a global flavour symmetry. 

The paper is organised as follows: In section 2, we review the holographic setup describing 
M = 2 SYM on an S 3 , studied in refs. [T3ll4H15j . We define a set of coordinates convenient 
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for the study of probe branes and their classification. Furthermore, we comment on the 
AdS/ CFT dictionary and the holographic renormalization of the probe branes, introduce 
the scaling exponents characterizing the self-similar behaviour of the theory near the topol- 
ogy changing transition of the D7-brane embeddings and review the calculation of the order 
of the phase transition [TJ]. 

In section 3.1 and 3.2, an external magnetic field along two of the S 3 directions is introduced. 
We construct the D7-brane embeddings and study the self-similar regime of the theory, ex- 
tracting complex scaling exponents. The analysis of scaling exponents is supplemented with 
numerical results to show that in an external magnetic field the phase transition is of first 
order. Next, we study the dependence of the fundamental condensate on the bare mass 
parameter for various magnetic fields. The effect of the magnetic field is to decrease the 
critical mass at which the confinement /deconfinement phase transition takes place. For suf- 
ficiently large magnetic field the critical mass vanishes and the transition happens at zero 
bare mass. Beyond this point there is no phase transition at all and at vanishing bare mass 
the stable vacuum features a non-vanishing negative condensate that spontaneously breaks 
the axial U(l) R-symmetry. To illustrate this, we construct a phase diagram of the theory 
summarizing this behavior. 

In section 3.3, the meson spectra of the theory are analysed. The effect of the magnetic field 
is to couple the vector and scalar modes. We show that the self-similar behaviour of the the- 
ory at finite magnetic field has a tachyonic instability. Across the phase transition, the meson 
spectrum displays a finite jump between the confined and deconfined phases. Moreover, we 
demonstrate that the globally symmetric vacuum develops a tachyonic instability for strong 
magnetic fields. However, it continues to be metastable for a small window of values of the 
magnetic field, even though the stable phase of the theory is the one with spontaneously 
broken global symmetry and negative condensate. The spectrum of the mixed modes ex- 
hibits Zeeman splitting at large bare masses, which leads to level crossing in the confined 
phase of the theory. For strong magnetic fields, the phase transition seizes to exist and 
the theory is in the confined phase. The ground state of the spectrum possesses a massless 
mode corresponding to the Goldstone boson associated with the spontaneous breaking of the 
global U(l) .R-symmetry of the theory, analogous to the rf meson in QCD. Furthermore, by 
studying the dependence of the mass of the meson on the bare quark mass near the origin, 
we find a characteristic M oc yfm behaviour reminiscient of the Gell-Mann-Oakes-Renner 
relation |18j . corroborating the existence of a massless Goldstone boson. 

2. Review of the model 
2.1. Flavours on S 3 

In this section we review the holographic setup describing flavoured M = 2 SYM on an S 3 , 
studied in refs. [T3~|14II15| . Let us consider the metric of AdS$ x S 5 in global coordinates: 

ds * = -(1 + r 2 /R 2 )dr 2 + r 2 dQ^ 2 + 1 + ^ + R 2 dQ 2 5 , (2.1) 
dtil = dej + cos2 M^s 2) 2 + sin 2 3 d^ • 
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It is convenient to define a new radial coordinate, 



u = ^(r + VR 2 + r 2 ) . (2.2) 



Then the metric (2.1) becomes 



Note that in this radial coordinate the metric has a conformal 1R 6 piece. However, u is 
bounded from below by u > R/2 and the transverse 1R 6 has a ball of radius R/2 at the 
origin. 

Next, we consider a stack of D7-branes embedded along the AdS 5 part of the geometry 
and wrapping an S 3 C S 5 . The dual field theory is Af = 4 SYM on S 3 coupled to Af = 2 
hypermultiplets. Note that, unlike the theory on flat space, a non- vanishing mass of the 
hypermultiplets breaks all supersymmetries and the theory preserves Af = 2 supersymmetry 
only for massless flavours. In the gravity setup this can be seen by studying the K-symmetry 
condition for the probe branes. 

The shape of the D7-brane embeddings can be determined by extremising the Dirac- 
Born-Infeld action, 

Sdbi = -Nffir J d 8 £e-*[-det(G ab + B ab + (27ra')F ab ] 1 / 2 . (2.4) 

Ms 

The radial part of the corresponding DBI lagrangian is given by, 

C K i 1 ~ m) i 1 ~ 5) 2 U " C ° S3 93 v/l+^W • (2-5) 

The resulting D7-brane embeddings fall into two classes [33]. Embeddings, which reach the 
origin of the AdS terminate on the shrinking S* 3 C AdSs. We will call these embeddings 
"ball" embeddings, because they end at u = R/2 in terms of the radial coordinate u. The 
second class of embeddings terminates above the origin of the AdS and instead wraps a 
shrinking S 3 G S 5 . We will refer to this class as "Minkowski" embeddings (cf. figure [l). The 
two classes are separated by a critical embedding which has a conical singularity located on 
the ball (the origin of AdS). The AdS/CFT dictionary relates the asymptotic behaviour of 
the transverse scalar 6^{u) to the vev and source of the fundamental bilinear, with the source 
being the bare mass of the hypermultiplet. The dictionary has been derived in ref. [19] using 
an appropriate holographic renormalization procedure. For large u one obtains the following 
expansion, 

m C\ m 

sin 3 = h -r - — ~ log u + . . . . (2.6) 

u u 2u A 

The condensate of the theory is then given by, 

(M) oc -2cj + m \og(m/R) = -2c , (2.7) 
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Figure 1: Blue curves correspond to "Minkowski" embeddings and Red curves correspond to "ball" 
embeddings. The black dashed line corresponds to the critical embedding. 



where we have defined a new parameter c proportional to the condensate. Note also that the 
bare mass of the hypermultiplet is proportional to m: the exact relation is m q = m/2ita' . 

By solving numerically for the D7-brane embeddings and studying the asymptotic be- 
haviour of 3 (u), one can construct a plot of the equation of state (i.e., a plot of c versus m). 
We generated this plot in figure |2j where we have used dimensionless parameters rh = m/R 
and c = c/R 3 . Note that there is no apparent multi-valued region near the transition from 
"Minkowski" to "ball" embeddings. In fact, by calculating appropriate critical exponents the 
authors of ref. [H] have shown that it is a third order phase transition. Furthermore, the au- 
thors argued that the phase transition corresponds to a quantum confinement /deconfinement 
phase transition triggered by the Casimir energy of the S 3 , which destroys the bound states 
of the quarks for sufficiently large curvature of S 3 . The state for which bound states exist 
corresponds to values of the parameter rh greater than the critical mass m* (the blue branch 
of the curve in figure |2|. The dissociated phase corresponds to values of the parameter rh 
smaller than m* (the red branch in figure [2]). The physical meaning of the parameter m 
is not obvious from its definition since the scale R represents the radius of S 3 , but is also 
related to the t'Hooft coupling A of the theory. The right way of thinking about m is as 
follows p^][T5] : rh = m/R = mR^/R 2 , where R3 denotes the radius of S 3 and R 2 oc y/X. 



(2.8) 



The exact expression for rh in terms of field theory quantities is 

7T m q R 3 



m 



V2 \/X 

For fixed bare quark mass m q , small values of rh correspond to small radius of S 3 and large 
Casimir energy triggering the dissociation of the bound quarks. Conversely, for large rh the 
radius of S 3 is large and the Casimir energy is small relative to the bare quark mass and 
thus the bound state is stable. 
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Figure 2: A plot of the condensate — c versus the bare mass fh. The states corresponding to 
"ball" embeddings are represented by the red curve and the states corresponding to "Minkowski" 
embeddings are represented by the blue curve. 



2.2. Critical exponents and the order of the phase transition 

In this subsection we focus on the vicinity of the critical embedding separating the "Minkowski" 
and "ball" classes of embeddings. We obtain the corresponding critical exponents and show 
(along the lines of ref. [14]) that the topology changing transition corresponds to a third 
order phase transition in the dual gauge theory. These studies will be applied in section 3 
to study the critical exponent of the transition in the presence of external magnetic field. 

To begin with, we zoom into the geometry near the tip of the critical embedding (the 
dashed curve in figure [T]). Let us consider the change of variables, 

u = 1/2(1 + z); y = n/2-e 3 ; (2.9) 



and expand the metric (2.3). To leading order in z and y, we obtain 

ds 2 zoom = -dr 2 + z 2 dQ<£ )2 + dz 2 + dy 2 + y 2 dQ (2) 2 + d<& , (2.10) 



which is just a flat metric on IR 1,9 . The Lagrangian for the D7-brane embedding (2.5) 
becomes 

3„3„ f\ I ni ,2 



c zoom = zYV^ + y' 2 ■ (2-n) 



The Lagrangian (2.11) is a special case of the general Lagrangian 



z 1 y 



k/2 T n ^l + y' 2 , (2.12) 



considered in ref. [20], for k = 6,n = 3. The solutions to the corresponding equation of 
motion have a scaling property y(z) — > ^y(fiz) and a critical solution y*(z) = \j2n/kz. 
Expanding in the vicinity of the critical solution y(z) = y*(z) + !;(z) results in the following 
equation of motion [20], 

z 2 e\z) + (n + k/2){z?(z) + £(z)) = . (2.13) 
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Equation (2.13) has a general solution 



£(z) = consti/ 1 + const 2 / 2 , (2.14) 

where 

5 1/2 = l/2 -(k/2 + n - 1) ± vW2 + n - l) 2 - 4(fc/2 + n) . (2.15) 

It is easy to see that for -2^ + 3 — n < k/2 < 2-\/2 + 3 — n the critical exponents <5i/2 are 
complex. In this window the solution exhibits a discrete self-similar behavior which seeds 
multi-valuedness in the equation of state and suggests that the phase transition is of first 
order (for more details, cf. refs [T4l2U|l21|l22|l23j ). In section 3 we will see that this holds 
true in the presence of a magnetic field. Interestingly without magnetic field we have real 
exponents [14"] . 

<Ji = -2; 5 2 = -3; (2.16) 
and the phase transition is continuous. Furthermore one can show that it is a third order 



phase transition [Ml . Note also that the linearized equation (2.13) is valid for large enough 



z (as the negative sign of the exponents suggests), while the zoomed in geometry (2.10) is 



valid for small z. Therefore our analysis is valid in an intermediate region of values for z and 



y, which exists since one can always consider sufficiently small constants in equation (2.14). 
Note that the D7-brane embeddings are uniquely determined by specifying initial conditions 
y(z ) = for "Minkowski" embeddings and z(y ) = for "ball" embeddings. We will focus 
on the "Minkowski" class of embeddings since the distance above the ball, specified by z , 
can be interpreted as a dynamical mass of the fundamental fields [2T] and acts as a natural 
order parameter for our phase transition. Now the scaling property of the equation of motion 
mentioned earlier suggests that if we rescale the initial condition by z' = Zq/ [a, the constants 



in equation (2.14) scale as: 

const^ = const i// 1-1 ; const 2 = const 2 // 2-1 ; . (2-17) 

In the field theory the critical embedding y* corresponds to a critical state characterized by 
some bare mass m* and fundamental condensate c*. Our next step is to assume that the 
field theory parameters m, c corresponding to the embeddings in the vicinity of the critical 
embedding depend analytically on the constants consti and const2- This is a reasonable 



assumption since the region where the linearized equation (2.14) holds is away from the 



conical singularity of the critical embedding. Therefore we can expand 

m — m, = A\ const! + A 2 const2 + .. . , (2-18) 
c — c* = B\ consti + B 2 const2 + . . . , 

where we have introduced new sets of constants A\, A 2 , B\, B 2 . Now using the scaling prop- 



erty (2.17) we can write: 



m'-m* = A 1 constx + A 2 const 2 /Z 2 " 1 + ... , (2.19) 
c' — c* = Bi consti /Z 1-1 + B 2 const2 // 2-1 + . . . . 
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Next, we eliminate the scaling parameter fi from equation (2.20) and obtain, to leading order, 



c - c* = DAm - raj + -D 2 (ra - raj 5 !- 1 + . . 



(2.20) 



where we have omitted the prime superscript and defined another set of constants Di,D 2 . 
In our case 5\ = —2, <5 2 = —3 and we find [H]: 



DAm — m*) + D 2 (ra 



ra* 



,4/3 



(2.21) 



Clearly, equation (2.21) suggests that the fundamental condensate of the theory is a contin- 
uous function of the bare mass with regular first derivative and divergent second derivative. 
Given that the condensate is a first derivative of the free energy, we conclude that the phase 
transition is of third order 



3. External magnetic field 

In this section we study the influence of an external magnetic field on the flavoured gauge 
theory. First we focus on the effect of the magnetic field on the confinement / deconfinement 
phase transition. For sufficiently strong magnetic fields we study the spontaneous breaking 
of the axial £7(1) -R-symmetry. We supplement our studies with an analysis of the meson 
spectrum. 



3.1. Introducing magnetic field 

In order to couple the fundamental fields to an external magnetic field we turn on a pure 
gauge 5-field along two of the directions of the S 3 where the dual gauge theory lives. Unlike 
the flat case considered in ref. [9] a constant S-field is not the natural choice on S 3 . It is 
instructive to write the metric of S 3 in terms of local tetrads: 

^ = e (l) 2 + e (2) 2 + e (3)2 (31) 

where the tetrads are defined by: 

e W =RdB x , e (2) = R sin 6 1 dfa , e (3) = R cos X dip! . (3.2) 
A natural choice for the -B-field is: 

B = He {l) A e {2) . (3.3) 

It is easy to check that the pure gauge condition B = dA is satisfied for A = —HR 2 cos a d/3. 
It is now straightforward to show that the DBI lagrangian from equation (2.5) is modified 
to: 

r 4 \ rrz # 2x : 




£ocwcos 3 # 3 ( 1-7^1 ) l-^a ) + H 2 R^l + u 2 6' 3 (u) 2 . (3.4) 

How does the external magnetic field affect the phase transition described in section 1? 
To answer that question we need to solve numerically for the D7-brane embeddings and 
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study the corresponding equation of state. However, at least for small magnetic field one 
would expect that the classification of the embeddings remains the same and there is still 
a critical embedding with a conical singularity at the origin of the AdS§ space (represented 
by the shell of the ball at u — R/2) separating the two classes of embeddings. It is then 
natural to calculate the corresponding scaling exponents. After performing the change of 



coordinates (2.9) and zooming into the vicinity of the critical embedding, we obtain the 



following Lagrangian 



C zoom oc R 2 Hzy 3 ^ + y' 2 ■ 
Notice that the power of z has changed and now we have the case k 



n 



(3.5) 

3 in equation 



(2.12). The critical exponents are then given by equation (2.15) and we obtain 



3 , .y/7 
2 2 



(3.6) 



Remarkably the critical exponents are complex^ Therefore the system exhibits discrete 
self-similar behaviour and the equation of state in the c versus m plane has a multi-valued 
behavior near the critical state (m*,c*), seeded by a spiral structure. Thus we expect that 
when we move away from the H = point, the third order phase transition becomes a first 
order phase transition. 

Following the same logic as for the continuous case considered in section 2.2, we consider 
approaching the critical embedding by scaling the initial conditions, namely z' Q = zo/fi for 
"ball" embeddings and y' = yo/p. for "Minkowski" embeddings. Then the corresponding 
values of the parameters rh and c scale as [23|20] : 



m' — m* 
c'-c* 



5/2 



M 



cos 



sin 



In jj, 



^-hip. 



sin 



cos 



^hi/i' 



2 



lnu 



M" 1 



m — m 
c- c* 



(3.7) 



where M is a constant non-singular 2x2 matrix. Equation (3.7) can be checked numerically 



by solving (numerically) the equations of motions derived form equation (3.4). In fact, it 
turns out to be more convenient to change variables, i.e., 



p = u cos #3 , L = u sin 9 3 



(3.8) 



The Lagrangian (3.4) becomes 



C oc p 



16(p 2 + L 2 ) 2 




1 - 



+ 



H 2 



4(p 2 + L 2 )J (p 2 + L 2 ) 



1 + 1/ 2 , 



(3.9) 



where we used the dimensionless coordinates L = L/R and p = p/R. The B-held introduces 
a new logarithmic divergence that can be regulated by adding an additional counterterm 
to the boundary, C n oc _B 2 /41og(p ma:c /-R). Note that his counterterm is independent of the 

fact these are the same scaling exponents as in the finite temperature case when the singular shell 
separating the two classes of embeddings is an event horizon. 
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bare mass and therefore does not modify the fundamental condensate of the theory. The 
parameters m and c are obtained by expanding L(p) for large p, 



r ~ ci m . 
L = m + — - — \ogp + ... 
p 2 2p 2 



(3.10) 



The condensate is calculated using equation (2.7), where the dimensionful parameters m, c 
and c\ are used. 

We now proceed to study numerically the self-similar structure of the theory near the crit- 
ical state (m*, c*) by approaching from the "Minkowski" class of embeddings. "Minkowski" 
embeddings are uniquely determined by specifying the initial value of L, namely, Lq = 1/(0). 
The critical embedding corresponds to Lq — L* = 1/2. We can consider some initial value 
Lq close to and then scale L' = Lq/ p. In this way we can trade the scaling parameter 



p in equation ( 3/71) for Lq/L' q and study m and c as functions of Lq for fixed m',c' and L' . 



This suggests that if we plot the quantities (m — m*)/(Lo — I/*) 5 ^ 2 and (c — c*)/ (Z — I/*) 5//2 
as functions of (a/7 /47r) log(L — L*), we should obtain trigonometric functions of unit pe- 
riod. In figure [3] we have presented our numerical results for H = 0.3. The dashed curves 
represent fits with trigonometric functions of unit period. One can see the excellent agree- 
ment as one explores states closer to the critical one (shifting toward negative values of 
(\/7/4f)log(L -L t )). 



[Lo-U 



2.5 . 2.0 J-L5 ' 1.0 • (» 



<1 
An 



log(io-L.) 



-2.0 t -1.5 .-1.0 



'. L5 

\'1.0 

\0.'5. 



-0.5 



-0,5 
-lip 
-1.5 



4^ 



log(Lo-L.) 



Figure 3: Plots representing the discrete self-similar structure of the theory near the critical state 
(m*,c*). The states closer to the critical one are to the left of the horizontal axis. The dashed 
curves represent fits with trigonometric functions of unit period. 



3.2. Equation of state and phase diagram 

In this subsection we analyze the dependence of the fundamental condensate on the bare 
mass. To this end we numerically solve the equation of motion for L(p) obtained from 
equation ( |3.9[ ) and extract the parameters m and c from the asymptotics of the solution at 
large p. 

In figure [4] we present our results for four different values of the magnetic field. As 
one can see, the effect of the magnetic field is to decrease the critical mass at which the 
confinement/deconfinement phase transition takes place. The dashed fitting curve in the 



10 




Figure 4: Plots of the condensate as a function of the bare mass for different values of the magnetic 
field. The black dashed curves represent a 1/rh fit. At H = H CT ~ 3.78 the critical mass vanishes. 
Beyond this point the theory is in a confined phase with broken chiral symmetry and the zero bare 
mass state has a non-zero negative condensate. 



figure represents a 1/rh fit for large rh. Moreover, for sufficiently large magnetic field H = 
H cr rs 3.78 the critical mass vanishes and the transition happens at zero bare mass. Beyond 
this point the phase transition seizes to exist and at vanishing bare mass the stable vacuum 
has non-vanishing negative condensate spontaneously breaking the axial U(l) R-symmetry. 
In our setup this breaking is analogous to the chiral symmetry breaking in the QCD vacuum. 
Thus we interpret this result as a manifestation of magnetic catalysis of chiral symmetry 
breaking in the flavoured gauge theory on S 3 . 

The shaded regions in figure [4] demonstrate the equal-area law which can be used to 
determine the critical mass. The condensate is a derivative of the free energy with respect to 
the bare mass and in this setup one can show that, by integrating numerically the condensate 
as a function of the bare mass, one can obtain the free energy (up to a additive constant) 
calculated by regularizing the euclidean on-shell action. This is why the use of the equal 
area law is justified. Note that from the plots in figure [4] one can observe that the disrete 
self-similar regime of the theory, analyzed in the previous section, is thermodynamically 
unstable. Our studies of the meson spectrum confirm that it is also unstable under quantum 
fluctuations. 

We proceed by obtaining the phase diagram of the theory. To this end we numerically 
generate the fundamental condensate as a function of the bare mass parameter for different 
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Figure 5: A plot of the phase diagram of the theory. For finite magnetic fields the phase transition is 
of first order, while at vanishing magnetic field the phase transition is of third order. For sufficiently 
large magnetic fields the phase transition seizes to exist and the theory is in a confined phase. 



values of the magnetic field. Then we integrate numerically to obtain the free energy and 
to determine the critical mass of the transition. Our findings are presented in figure [5] The 
phase diagram is similar to the one studied in ref. [2~jlT2"5] . The analogue of the temperature 
in our scenario is the Casimir energy of the theory on S 3 , which dissociates the mesons 
states. In our case however, the phase transition takes place at zero temperature and is thus 
a quantum phase transition. An interesting property of the phase diagram is that the first 
order phase transition for finite magnetic fields becomes a third order phase transition for 
vanishing magnetic field. One can also see that for sufficiently large magnetic field there is 
no phase transition at all and the theory is in a confined phase. 

3.3. Meson spectra with external magnetic field 

Another focus of the present work is to study the effect of an external magnetic field on 
the spectra of mesons in global AdS. The presence of a magnetic field manifests itself, as 
expected, by Zeeman splitting and level crossing. This will be studied in detail below. We 
explore the dependence of the spectra on the bare mass parameter for a large range of values 
of the magnetic field. 

3.3.1. Derivation of the fluctuation equations 

To investigate the spectrum of light mesons, we study quadratic fluctuations of the D7-brane 
embedding along the transverse coordinates L, 03, which we expand in the following way: 

L(p) = L (p) + 2W X (p), MP) = + 2W$(p), (3.11) 
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where Lq(p) is the classical D7-brane embedding. The discussion closely follows earlier work 
reported in [9.24 25 26J. Recall the induced metric of the D7-branes in (L, p) coordinates, 



ab D7 



+ 



p 2 + L 2 
R 2 
R 2 



1 + 



R 2 



4:(P 2 + L 2 ) 



2 , P 2 + L 2 



dr 2 + 



R 2 



R 2 



4(p 2 + L 2 ) 



(1)2 



p 2 + L 2 



(i + L'( P ) 2 )dp 2 + p 2 dn 



(2)2 



The equations of motion can be obtained straightforwardly from the appropriate DBI and 
WZ actions: 



S — Sbbi + Swz, 
SvBi = -N f p 7 J d 8 £e-*[-det(E ab + 2ira'F ab )} 



1/2 



Ms 



Swz = I F {2) A F {2) A P[C (4) ] + (2na>)p 7 



F{2) A -B( 2 ) A P[C(4)], 



(3.12) 



where P [C^] is the pull-back of the 4-form potential sourced by the stack of N c D3-branes, 
and P[C( 4 )] is the pull-back of its magnetic dual. Expanding E ab to second order in a', we 
have 

E ab = E° ab + + (2rra') 2 E 2 ab , (3.13) 

where 



E ab 



9ab(L (p)) + B ab , 
El b = G LL L' (p) (d aX 5 p b + d bX 5P) + (d Lo9ab ) X , 

E 2 ab = G LL {d aX d bX + L 2 d a ^) + (d Lo G LL L' ) (d aX 6 p b + d bX 5 p a ) X + \d 2 u G abX 2 , 

where g ab and B ab are the induced metric and B-field on the D7-brane world volume and the 
prime stands for derivation w.r.t. p. Splitting (E® b ) 1 into symmetric and anti-symmetric 
parts, i.e., (P^)" 1 =: S ab + J ab , where 



S ab = diag i -g a 



933 



933 



-1 „-l „-l 



-1 



J' 



ah 



H 



9h + H 2 



g2 _|_ jj2 > g2 _|_ jj2 > ^33 ' 9 PP ) 9e 2 e 2 ' 9<f, 2 fa > 9^ 2l p 2 ( > 
{5 a 2 5 b - S\6l) , 



a - P 2 + L »(l 
9u - [I 



R 2 



4(P 2 + LI) 



933 



P 2 + L 



R 2 



1 - 



R 2 



4(P 2 + LI) 



9pp 



R 2 



P 2 + Ll 



(1 + L 



It was demonstrated in refs. PI2"T] that the effect of the magnetic field regarding the equations 
of motion is to couple the scalar and vector modes; namely, $ will couple to the A and A 3 
components of the gauge field, while the X fluctuations will mix with the A\ and A 2 in the 
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presence of a magnetic field. This can be understood from the fact that our ansatz for the 
external magnetic field breaks part of the remaining 3+1 dimensional symmetry. Note that 
Lorentz invariance is broken already, cf. (2.3). 

With this in mind, we arrive at the following expressions for the relevant terms in the 
Lagrangian: 

S tt 

—<), \ a, \ ■ ■■ ^ , 



2 9{p)Gll 1 + L 



1 S pp 
,2 d tXdtX+ ^g{p)G L L ^ , r/2 d p xdpX 



1 + ^ 2 



(3.14) 



+ 



L' 



d Lo (d Lo log^(p)) - - + ° L ,2 d p ( d L ^gg(p)) 



X 



(2) 
AA 



Cf A = -H{d p K{p))$F { 



03 • 



The metric component G<p 3r f, 3 and the functions g(p) (the Lagrangian density) and K(p) are 
given by 



G 



<p3<t>3 



R 2 Ll 



P 2 + LV 
g(p) := yf-detEft, = p 3 (l 

K{p) 



R 4 



16(p 2 + L 2 ) 2 



R 2 



4(p2 + L \ 



H 2 R 4 



4 „4 



R 4 p 



(P 2 + L 2 ) 2 ' 



1 + ^ (P) 2 , 



(3.15) 



Now it is straightforward to obtain the equations of motion for the fluctuations from the 
usual Euler-Lagrange procedure, yielding 



9{ P ) dp 



g(p)d P x ' 



R l ^x 



(1 + L' 2 ){p 2 + LI) 2 (l + 



R 2 



d Lo (d Lo log g(p)) - 2 d p (d Lo log g(p)) 



m 9 ' 



W 9 ' 



'g(p)Lld p $ 



TI2 

^0 



1 + L' C 
L 2 R 4 u 2 $ 



(P 2 + L 2 ) 2 (l + 



Hd p K 

-re 



R 2 



*(P 2 +L 2 ) 



g(p) 



x = o, 
o. 



03 



g(p)d P F 03 



R 2 



+ 



R 4 



(P 2 + L 2 ) 2 1 



16(p 2 +Lg) 2 



cu 2 F 03 - 



Hd p K 



(3.16) 
(3.17) 

V$ = o, 



(3.18) 



where F 3 = doA 3 — d 3 A . We have assumed a plane wave ansatz for the fluctuations of 
the form SX(t,p) = e~ tult 5X(p). A few remarks are in order: We will be interested in in- 
vestigating the spectrum of "pions" which correspond to fluctuations along 3 . Therefore 
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we will only consider the gauge field components Aq and A% which couple to $ and solve 
the corresponding coupled system of differential equations. Similar equations of motion were 
obtained and discussed e.g. in p3] for the case of AdS 5 x S 5 . In order to be able to derive 
an equation in terms of F 3, we have to define F 03 = d T A 3 , setting d 3 A = 0. This is due to 
the broken 50(1, 1) symmetry of the Minkowski part of (2.3). Moreover, since we assume 
diX = 0, the x equations decouple from those for the gauge field components A\ and A 2 . 
In the following sections, we will present numerical solutions to the equations of motion 
(3.16 )-(3.18 ) which were obtained employing a shooting technique in Mathematica. The 
general strategy is to start with appropriate initial conditions and shoot towards the bound- 
ary where we are interested in finding those solutions that display the expected fall-off 
behaviour in the UV. 



3.3.2. Fluctuations along L 

Here, we will present our numerical study of the meson spectrum associated with fluctuations 
along L. As a first step, we obtain the spectrum at zero bare mass, i.e., for the trivial 
embedding wrapping the S 3 within the S 5 . 

Zero bare mass. Figure [6] shows the dependence of Ui on the strength of the magnetic 
field H for the lightest four meson states, i — 1, ... ,4.. We plot sgn(o; J 2 )|a;j| which becomes 
negative when Ui becomes imaginary (i.e., the state becomes tachyonic). For H — > 0, the 

sign[tl> 2 ] |<31>| 




Figure 6: The spectrum of x fluctuations for zero bare mass (trivial embedding) vs. the magnetic 
field H. 



spectrum is discrete and equidistant with eigenfrequencies given by (15] 

u = (2n + 3)^-. (3.19) 
R 

Increasing the magnetic field, we observe that the lightest meson state becomes tachyonic, 
i.e., u\ becomes imaginary, for some H CI > H*. This indicates the existence of a metastable 
phase between H* and H CI , which becomes unstable for even larger H. Note that only the 
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scalar modes condense, i.e become unstable, whereas the vector modes remain stable. This 
is not in contradiction with the proposed condensation of the charged p meson discussed in 
|53j . since we consider only one type of flavour and therefore all meson modes are neutral. 
Nonzero quark mass. We proceed by investigating the spectrum of fluctuations along 
L as a function of the quark mass parameter fa. We again solve numerically the equation 
of motion for Xi ec L- (3.16), for the two classes of embeddings, "Minkowski" and "ball" 
embeddings. 

For intermediate values of the magnetic field we present our results in figure [7} The 
dashed lines in the figure represent the spectrum of the theory on R 1 ' 3 without external 
magnetic field, studied in ref. [28J, 



2m 



w = ^-V(n + l)(n + 2), 



(3.20) 



As expected, at bare masses larger than the energy scales set but the magnetic field and 



the Casimir energy, the spectrum is well described by equation (3.20). As one approaches 



H=0.1 



H=2.0 




Figure 7: The spectrum of \ fluctuations vs. bare quark mass rh for intermediate values of H, 
H = 0.1 for the left plot and H = 2.0 for the right one. The vertical dashed line indicates the 
location of the critical bare mass m* , while the coloured dashed lines correspond to the spectrum of 
pure AdSs x S 5 . The red curves correspond to the ball embeddings and the blue curves represent 
the Minkowski embeddings; the phase transition happens for the critical embedding separating the 
two classes. 



the phase transition, the spectrum becomes multivalued, with "competing" confined and 
deconfined phases. At the phase transition the spectrum has a finite jump. It is also 
interesting to see that close to the critical state, where the theory has a discrete self-similar 
structure, the spectrum becomes tachyonic, which suggests that the self-similar regime is 
unstable under quantum fluctuations and cannot be realized by "supercooling" (i.e., it is not 
meta-stable). In the deconfined space the spectrum remains discrete because the theory is 
in a box. 

For magnetic fields above criticality, (H > H* « 3.78), the only stable phase is the 
confined phase. The corresponding spectrum is presented in figure |8j As one can see the 
positive m branch is stable. 
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H=4.0 

sign[<y 2 ]|a>| 




0.5 1.0 1.5 2.0 2.5 3.0 3.5 



Figure 8: A plot of the spectrum of fluctuations along L for H = 4.0 > H*. One can see that the 
positive rh branch is stable. 



3.3.3. Fluctuations along fa 

The spectrum of scalar mesons associated with fluctuations along fa is particularly interest- 
ing because it possesses a mode analogous to the r/-meson of large N c QCD. The analogy 
comes form the fact that rotational symmetry along 03, present when L = 0, corresponds to 
the anomalous axial U(l) R-symmetry, which is restored at large N c . Spontaneous breaking 
of this symmetry leads to a massless Goldstone boson, the analogue of rf in QCD. In the 
gravity setup the repulsive potential due to the S-field decreases the asymptotic separation 
L(oo) of the D7-brane embeddings, and there exists an embedding which asymptotes to 
zero separation at infinity but has finite separation in the bulk of the geometry. Thus, the 
rotational symmetry along 03 present at infinity (L(oo) = 0) is broken in the bulk, which 
corresponds to a spontaneous breaking of the dual U(l) .R-symmetry [29]. By generating 
the spectrum of fluctuations for various values of the bare mass parameter rh we will show 
that indeed the ground state is massless at vanishing bare mass. Furthermore, for small 
values of rh we will demonstrate a characteristic Gell-Mann-Oaks-Renner \frh q dependence 
[1~8] of the spectrum. 

In the analysis, we have to take into account the non-trivial mixing of the <J>-mode with 



the gauge field components A ,A 3 , or equivalently F03, cf. (3.17). Again, we will study the 
meson spectrum numerically and require normalizability of the solutions in the UV. This 
condition will again lead to a discrete spectrum. 

The numerical results are presented in figure [9j where we plot u versus rh. For large 
bare quark mass rh, the spectrum is expected to match the spectrum of the flat D3/D7 



intersection on AdS^ x S 5 , given by equation (3.20), and our results in the asymptotic region 



are indeed in good agreement with that expectation. As before, the red and blue curves 
correspond to the "ball" and "Minkowski" embeddings, respectively. One can see that near 
the phase transition (represented by the vertical dashed line in figure P), the spectrum is 
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H=0.35 



15 



10 





0.5 i 1.0 1.5 2.0 

Figure 9: The spectrum of $ fluctuations vs. quark mass fh for various values of H. The 
vertical dashed line indicates the location of the critical bare mass m*, while the coloured 
dashed lines correspond to the spectrum of pure AdS^ x S 5 . The red curves correspond to 
the ball embeddings and the blue curves represent the Minkowski embeddings; the phase 
transition happens for the critical embedding separating the two classes. 



multivalued and at the phase transition it has a finite jump. However, unlike the spectrum of 
fluctuations along L, there are no tachyonic modes in the self-similar region near the critical 
state. As expected, our investigation confirms a Zeeman-like effect, namely a splitting of 
states (which will be proportional to the magnitude of H at weak magnetic field or large fh). 
This can be gleaned from figure [9] where we observe two separate lines emanating from each 
asymptotic meson state at large bare mass fh. At smaller values of fh the Zeeman splitting 
is strong and the energy levels intersect. This phenomenon is known as "level crossing". 



Gell-Mann-Oakes-Renner (GMOR) relation. In figure 10 we have presented the spec- 
trum of fluctuations for strong magnetic field (H = 4.0 > if*) when the phase transition 
disappears and theory is always in a confined phase. The spectrum is again featuring Zee- 
man splitting and intersection of energy levels (level crossing). At vanishing bare mass the 
spectrum has a massless Goldstone mode. Zooming into the small bare mass region, the 



bare mass dependence of the ground state meson mass (the right plot in figure 10) shows 



the characteristic Gell-Mann-Oaks-Renner relation, M oc vm, cf. [26] . 
4. Conclusions 

In this paper we studied the influence of an external magnetic field on a flavoured large N c 
gauge theory on S 3 . We find that there is a competition between the effect of dissociation 
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of the meson states due to the finite volume Casimir energy of the theory and the effect 
of the magnetic field which favours bound-states of mesons. As a result, the theory has 
an interesting phase diagram consisting of confined and deconfined phases separated by a 
critical curve across which there is a first order confinement /deconfinement quantum phase 
transition. At vanishing magnetic field the critical curve ends in a point where the phase 
transition is of third order [14.]. For sufficiently strong magnetic fields, the phase transition 
seizes to exist and the theory is in the confined phase. In this regime the vacuum spon- 
taneously breaks the global U(l) -R-symmetry by having a non-zero negative fundamental 
condensate. This is an example of magnetic catalysis of chiral symmetry breaking of a U(l) 
axial symmetry. Thus, the associated Goldstone boson is the analogue of rf in QCD. The 
effect of the external magnetic field on the meson spectra is to couple the scalar and vector 
modes. At large bare masses the spectrum exhibit Zeeman splitting of the energy levels, 
which leads to level crossing in the confined phase of the theory. Across the phase transition 
the spectrum has a finite jump between the confined and deconfined phases. For sufficiently 
strong magnetic fields the only stable phase is the confined phase and the ground state of 
the spectrum possesses a massless mode corresponding to the Goldstone boson associated 
to the broken global U(l) i?-symmetry of the theory, in analogy to the rf meson in QCD. 
Furthermore by studying the dependence of the mass of the meson on the bare quark mass 
near the origin we have found a characteristic M cx y/m Gell-Mann-Oakes-Renner relation. 
A possible extension of our results would be to study the scenario with magnetic field at 
finite temperature. However, at finite temperature the phase transition is thermal and the 
finite volume does not change the qualitative behaviour of the theory. Thus we expect the 
theory to have similar qualitative behaviour as in flat space. 

An interesting direction for future work is to consider the effects of a magnetic field and var- 
ious chemical potentials on fields theories with and without defects on compact manifolds. 
Defect field theories can be realized in an holographic framework, e.g., by the introduction 
of D5-brane probes [T4fl5] . 
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